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Abstract 

Through the hysteresis loop and magnetization spatial distribution we study and 
compare two models for surface anisotropy in nanomagnets: a model with transverse 
anisotropy axes and Neel's model. While surface anisotropy in the transverse model 
induces several jumps in the hysteresis loop because of the cluster- wise switching 
of spins, in the Neel model the jumps correspond to successive coherent partial 
rotations of the whole bunch of spins. These calculations, together with some hints 
from available experimental results, suggest that Neel's model for surface anisotropy 
is more appropriate. 
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1 Introduction 



Understanding surface effects in magnetic nanoparticles is a big challenge from 
both the experimental and theoretical points of view. Experimentally, this 
challenge resides mainly in the fact that the current experimental techniques 
do not allow for probing the microscopic structure of the surface in round 
geometries, as is the case for 2d magnetism [1]. In particular, we so far have 
no clear idea as to what kind of surface anisotropy we have in a nanoparticle. 
The corresponding constant has been estimated by many authors in an indirect 
way from neutron scattering experiments [2] , magnetic measurements [3] , and 
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FMR experiments [4], by fitting the results assuming an effective anisotropy 
constant. This is a rough estimation that is only valid when deviations from 
collinearity are very small. On the other hand, the easy direction of the magne- 
tization on the surface, let aside that of the atomic magnetic moments, is still 
unknown. In fact, even the anisotropy in the core of the particle is assumed 
to be that of the bulk material. As such, to investigate surface effects on the 
properties of nanoparticles theorists can only "try" some models. The most 
often used model is what we call here transverse surface anisotropy (TSA) 
model which attributes an easy axis along the transverse (radial) direction to 
each spin on the boundary [see the textbook [5] and references therein]. An- 
other more physically plausible model is what we call Neel surface anisotropy 
(NSA) model, developed by Neel [6]. This model is more realistic because the 
anisotropy at a given lattice site only occurs if the site's environment presents 
some defects, e.g., in the coordination number. In [7] (see also [8]) we used 
the TSA model and studied the details of the switching process in a spherical 
particle when the surface anisotropy constant K s is varied. We solved the cou- 
pled Landau-Lifshitz equations with local and global rotation constraints. In 
particular, we found that there is a characteristic value of K s separating two 
different regimes: For small K s , as compared to the exchange coupling, the 
switching is basically Stoner-Wohlfarth like, i.e., coherent, whereas for large 
values the magnetization switching can no longer be regarded as coherent, as 
it operates cluster-wise. In [9] we used the Green's function technique in the 
continuum limit, in the absence of magnetic field, and computed for the NSA 
model the contribution of the surface to the anisotropy energy of a spherical 
particle with simple cubic lattice structure. We found that this contribution 
is second order in the surface anisotropy constant K s [see Eq. (7) below], 
scales with the particle's volume M and has cubic symmetry with preferred 
directions [±, ±,±]. 

The present work follows on the previous ones [7] , [9] and is an attempt to give 
some hints for answering the question addressed in the title, namely which one 
of these models is most adequate to describe surface effects in a small mag- 
netic system. We have no pretention here to give a definite answer to this 
question as we believe that only with the help of experiments, e.g., measur- 
ing the hysteresis cycle of an isolated nanoparticle, could one do so. However, 
today technical difficulties prevent such measurements [see the review [10]]. 
Meanwhile, it is important to understand the effects of surface anisotropy 
within the framework of the current models, and eventually make some pre- 
dictions for future investigations. In addition, we do not try to be exhaustive 
since a detailed study of the magnetization switching within the TSA model 
has already been presented in [7] . Instead, we only investigate the main differ- 
ence between the two models with regard to the hysteresis loop and magnetic 
structure so as to gain some insights into their relevance to the study of surface 
effects in nanomagnets. More precisely, we compute the hysteresis loop and 
magnetic structures of a spherical particle with simple cubic lattice structure 
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and uniaxial anisotropy in the core with easy axis z. We consider both cases 
of vanishing and non-vanishing exchange interaction between spins inside the 
particle. The unrealistic non-interacting case is considered only because it 
allows for a comparison of the TSA and NSA in pure form. 



2 Statement of the problem 



Our model Hamiltonian is given by the (classical) anisotropic Dirac-Heisenberg 
model [7] 

N 

H = - J2 J ij S i ■ S j - (9Vb)H ■ J2 S i + H <™, (1) 

(ij) i=l 

where Sj is the unit spin vector on site i, H the uniform magnetic field, M 
the total number of spins (core and surface), and Jij{= J > 0) the nearest- 
neighbor exchange coupling. H an can be the uniaxial single-site anisotropy 
energy 

H an = -Y J K l {s l -e l )\ (2) 



with easy axis e^ and constant Ki > 0. If the spin at site i is in the core, the 
anisotropy axis e^ is taken along the reference z axis and Ki = K c . Otherwise, 
for surface spins, this axis is along the radial (or transverse) direction and 
Ki = K s . In this case, the model in (2) for surface spins is the TSA model. 
The second model that will be studied and compared with the previous one, is 
the more physically appealing model of surface anisotropy that was introduced 
by Neel [6], 

tf el = ^EEW 2 , (3) 

% j=i 



where Z{ is the coordination number of site % and = r^jr^ a unit vector 
connecting the site i to its nearest neighbors j. This model is more realistic 
because the anisotropy at a given site occurs only when the latter loses some 
of its neighbors, i.e., when it is located on the boundary. This is the NSA 
model. 

In this work we ignore dipolar interactions between spins inside the particle, 
since in our case of small spherical particles, the volume term is negligible [11] 
and the shape anisotropy yields an irrelevant constant. In addition, here we 
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restrict ourselves to the investigation of "pure" surface anisotropy within the 
two models. 

In the sequel, we will use the following notation h = H/2K C for field, with the 
magneto-crystalline ainsotropy constant K c in the core taken as the energy 
scale, j = J/K c for the exchange coupling, which will be taken in our calcu- 
lations the same everywhere inside the particle, k s = K s /K c , and D denotes 
the particle's diameter. 

We study hysteretic properties of a nanoparticle by solving the (local) Landau- 
Lifshitz equations at zero temperature. The method used to obtain the hys- 
teresis loop and other characteristics were given in great detail in [7] and will 
not be repeated here. We compare the implications of the two models, TSA 
and NSA, on the hysteretic properties of a nanoparticle, such as hysteresis 
loop and coercive field as a function of the surface anisotropy strength k s . In 
all these calculations the anisotropy in the core is taken as uniaxial along the 
reference axis z, and the field is applied at an angle of 7r/4 with respect to the 
latter in the x — z plane. 



3 Results and discussion 

3.1 Noninteracting case 

In order to compare the TSA and NSA models in pure form, we first consider 
the case of a spherical particle of non-interacting spins (Jy = 0), and the 
core-anisotropy constant equal to that on the surface. In Fig. 1 we plot the 
hysteresis loop for both TSA and NSA models. First, in both cases the core 
switches at the same field of 1/2 since the anisotropy axis in the core is along 
z and makes the angle 7r/4 with respect to the field direction. In the case of 
TSA, at zero field the surface spins are obviously directed along the radial 
direction [see Fig. 2 h — 0]. At non zero field, as was discussed in detail in [7], 
because of the radial direction of the surface anisotropy easy axes, the latter 
make different angles with respect to the field direction, and hence switch at 
different values of the applied field. This leads to a cluster-wise progressive 
switching of surface spins, as can be seen in the structures of Fig. 2 at the 
respective field values h = 0.47,0.56,0.64, 1.12 at which no spins, 2 spins, 3 
spins, and all spins have switched. For the surface anisotropy constant used 
in Fig. 2 the core switches before the full switching of the surface. 

In the case of NSA, the situation is quite different. The hysteresis loop exhibits 
a first jump at a very small field, and more importantly this jump corresponds 
to the switching of the entire surface. This can be seen in Fig. 3 by examining 
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Fig. 1. Hysteresis loops for a spherical particle with D = 10 (M = 360), j = 0, k s = 1 
for TSA and NSA models. 

the structures at the field values h — and ft, = 0.04. For smaller particles with 
more spins on the boundary than in the core, the first jump would also cor- 
respond to the switching of the particle's magnetization. Moreover, in Fig. 1 
there appear two jumps corresponding to surface spins with two different local 
environments: the spins on the faces and those located near the corners. How- 
ever, this difference in the switching field is just a numerical artefact, since as 
will be shown shortly these two groups of spins jump at the same field value, 
that is zero. Indeed, a spin on the face has 5 nearest neighbors, so that the 
local energy given by (3) yields, 



e^2§- = -(4 + 4 + ^) + ^- h - s * W 

— > - cos 2 9 -= [sin 6 cos y? + cos 6} . 

2 v 2 

where 9, tp are the spherical coordinates of the spin s. The arrow in Eq. (4) 
indicates that the irrelevant constant s 2 + s 2 + si = 1 is discarded. So, in zero 
field the energy (4) has a minimum in the x — y plane. Similarly, a spin near 
the corner has three nearest neighbors (two on the same plane and one on the 
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Fig. 2. Magnetic structures corresponding to the middle plane of the particle for the 
TSA model. The field values are indicated on top of the structures and correspond 
to the ascending hysteresis-loop branch of Fig. 1 for TSA. 

adjacent plane), so that its energy reads 



/ 2 i 2 i 2\ 

e * = -^Ks x + s y + s z ) 



h • Si 



_h_ 

72 



[sin 9 cos 4> + cos 9] 



(5) 



whose minimum, in zero field, is along the diagonal of the cube [see Eqs. (7, 
8) below and Ref. [9]]. 

As noted above, while the smallest switching field for a TSA spin is 1/2, here 
we see that the switching field of an NSA spin may assume very small values. 
Indeed, for spins on the faces and near the corners, it can be shown from 
(4) and (5), respectively, that the switching field is zero. For example, the 
magnetization direction which yields the lowest energy of an NSA spin on the 
face is given by if = and 



cos 9 n 



h_ 

V2' 



(6) 
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Fig. 3. Magnetic structures corresponding to the middle plane of the particle for the 
NSA model. The field values are indicated on top of the structures and correspond 
to the ascending hysteresis-loop branch of Fig. 1 for NSA. 



which when inserted in the determinant of the energy Hessian yields a zero 
switching field. This is confirmed by the result of the numerical calculations 
in Fig. 1, up to numerical precision. 



In fact, the reason for such small switching fields for surface spins in the NSA 
model is quite simple, and one actually does not need the above calculations 
to convince onself. Indeed, for spins on the face there is a rotation symmetry 
which makes all directions in, e.g., the x — y plane degenerate, see Eq. (4). So, 
when the field increases from the negative saturation value towards zero, these 
spins select the direction given by the projection of the field direction e h on 
the x — y plane. As soon as the field points into the opposite direction, even a 
very small value thereof is sufficient to make the spins on the face reverse their 
direction. The same arguments hold for spins near the corners, since these too 
have full symmetry in the plane containing their easy direction. 
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3.2 Interacting case 



In order to investigate the surface influence on the core and vice-versa, we 
consider the same particle as before but now with realistic non-vanishing ex- 
change coupling everywhere in the particle. The results of the hysteresis loop 
and some of the corresponding spatial distributions of the magnetization are 
shown in Figs. 4 and 5, respectively. In this case of fully interacting spins, we 
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Fig. 4. Hysteresis loops for a spherical particle with D = 10, j = 1, k s = 0.15. 

see that as the field goes from to 0.11, the magnetization of a NSA particle 
experiences its first jump, while that of the TSA particle only decreases (in 
absolute value) in a progressive manner. This is due to the fact that the angle 
of rotation between the direction at h — and that at h = 0.11, and thereby 
the difference in the projection of magnetization along the field, is much larger 
for the NSA particle than for the TSA one. The reason is that at zero field, in 
the case of the TSA model, all core spins and those surface spins located on 
the x — y faces are directed along the easy axis z, and hence make up a strong 
effective field that forces the spins on the x — z,y — z faces into nearly the 
same direction as in the core. As to the NSA model, there are near-corner spins 
directed along the diagonal of the cube, spins on the x — y faces which lie in 
these faces, and finally the more numerous core spins directed along their easy 
axis z. This results in a competition which ends up with a net magnetization 
in a direction inside the cube with its largest component along z. These results 
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Fig. 5. Magnetic structures for the TSA and NSA particles with 
D = 10, j = l,k s = 0.15. The field values on top of each structure correspond 
to the ascending hysteresis-loop branches of Fig. 4. 

suggest that the surface has a stronger effect in the case of NSA model as com- 
pared to the TSA model. Indeed, in Fig. 5 (NSA), for instance, we see that at 
zero field the core spins are deviated from their easy direction. Moreover, in 
the non-interacting case, in the NSA model, the surface jumps in a coherent 
way at a much lower field than the core which also switches coherently. In the 
interacting case, the surface imposes, via exchange coupling, its switching at 
low fields on the core, and this results in an abrupt partial coherent rotation 
of all spins, core and surface. The core then proceeds to switch and drives the 
whole bunch of spins at a higher field. 

In [9] we studied the NSA model using both the numerical solution of Landau- 
Lifshitz equation at zero temperature and the analytical Green's function tech- 
nique in the continuum limit, though preserving the discreteness of the lattice 
by introducing a spin density function. We showed that the contribution of 
the NSA to the energy of a spherical particle, in the absence of core anisotropy 
and applied field, reads, 



£ = K 



Jo 



+ m 4 z ), 



(7) 



where m a ,a = x,y,z, are the components of the net magnetization of the 
particle, Jo = z J is the Fourier transform of the exchange coupling J^, z = 6 
the coordination number in the core, and k a surface integral. Computing the 
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Fig. 6. Hysteresis loop for an NSA particle with j = 1, D = 10 and variable k s 
energy difference between two major orientations of m leads to 
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(8) 



This indeed confirms what we said above, namely that NSA has a stronger 
effect than TSA. 

In Fig. 6 we plot the hysteresis loop as obtained from the NSA model for 
different values of the surface anisotropy constant k s . These results demon- 
strate that actually the amplitude of the abrupt coherent rotation and the 
field at which it occurs increase with k s . Moreover, the length of the subse- 
quent "shoulder" corresponding to a smooth rotation of the whole spins and 
that precedes the magnetization switching decreases with k s and vanishes at 
some higher value of k s (0.75 here). This means that as k s becomes large, 
the amplitude (or angle) of the abrupt coherent rotation, induced by surface 
anisotropy, becomes so large that this rotation coincides with the switching of 
the net magnetization and thus no other steps are observed. 
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In summary, surface anisotropy in the TSA model induces several jumps in 
the hysteresis loop because of the cluster-wise switching of spins [7]. This 
is due to the single-site nature of the anisotropy in this model, and to the 
radial direction of surface anisotropy axes, and thence to the distribution of 
angles that these axes make with the field direction. In the NSA model, the 
existence of several jumps in the hysteresis loop does not require high values of 
the surface anisotropy constant [see Fig. 4], and is due to successive coherent 
partial rotations of all spins in the particle. For large values of the surface 
anisotropy constant, there is only one big coherent rotation which coincides 
with the net magnetization switching. 

As in Ref. [7] where the switching field was also studied as a function of the 
surface anisotropy constant k s , we have performed the same calculation here 
for the NSA model, for the sake of comaprison. The results are presented in 
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Fig. 7. Switching field against surface anisotropy for a spherical particle with 
D = 10,j = l,ip = 45, for both TSA and NSA models. The inset gives h c (k s /j) for 
NSA alone. 

Fig. 7. For the TSA model, there is a sudden increase in h c (k s /j) occurring 
at a value of k s /j of order one. This was interpreted in [7] as a critical value 
marking a departure from the Stoner-Wohlfarth model, i.e, a transition from 
the coherent-rotation regime, where the hysteresis loop exhibits only one jump, 
into a regime of cluster-wise magnetization switching, where several jumps 
appear in the hysteresis loop. In the case of NSA, h c assumes much lower 
values than in the TSA model as already discussed earlier, and increases in 
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a more progressive way. On the other hand, it is apparent from the inset in 
Fig. 7 that there is in fact some kind of transition between different regimes, 
though with a lower "critical value" for k s /j than for TSA model. 



4 Conclusion 

We have compared the effects of transverse surface anisotropy (TSA) and 
Neel's surface anisotropy (NSA) on the hysteretic properties of a spherical 
particle. In the non-interacting case, i.e., with zero exchange coupling every- 
where inside the particle, the NSA induces more jumps than the TSA in the 
hysteresis loop the first of which occurs at very low fields, and corresponds 
to the switching of the entire surface. This is due to the rotational symmetry 
which renders all directions in, e.g., the x — y plane, degenerate. Moreover, 
in the NSA model, the spin switching within the surface and within the core 
is coherent. In the interacting case with the same exchange coupling in the 
whole particle, we still have qualitatively the same hysteresis loop as before, 
but now the first jump in the NSA model occurs at a field slightly larger than 
before, and more importantly corresponds to a coherent partial rotation of all 
spins (core and surface). More precisely, in this case the surface imposes its 
switching on the core, and thereby the two jumps in the hysteresis loop, in 
fact, correspond to two successive coherent rotations of the whole bunch of 
spins inside the particle. This shows that, for the same constant, the surface 
anisotropy in the NSA model has a stronger effect than in the TSA model. 
Moreover, spin switching in TSA is cluster-wise, while that for NSA is coherent 
but operates by successive coherent partial rotations. 

Experimental results on individual cobalt clusters [see the review [10]] on the 
Stoner-Wohlfarth astroid, i.e., the angular dependence of the switching field, 
together with the switching field distribution and the non-switching proba- 
bility, have confirmed the validity of the Stoner-Wohlfarth and Neel-Brown 
models describing the magnetization reversal by uniform rotation [12] at zero 
and non-zero temperature, respectively. On the other hand, even though rela- 
tively large surface anisotropy constants are obtained for cobalt particles [13], 
the results are still well fitted by the aforementioned models that are based on 
coherent rotation of the net magnetization. This is at variance with what one 
expects from the TSA model with large anisotropy constant, as was demon- 
strated in [7] and reviewed here. Therefore, in addition to the fact the NSA 
model is intuitively more realistic, these experimental indications together 
with the present calculations give a further hint that the surface anisotropy in 
nanomagnets is most probably of the NSA kind. However, a more convincing 
check would only come from direct measurements of the hysteresis loop of an 
individual nanoparticle, were they to become possible. 
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